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ABSTRACT 

We present a measurement of the X-ray spectrum of the luminous X-ray binary in I Zw 18, the 
blue compact dwarf galaxy with the lowest known metallicity. We find the highest flux yet observed, 
corresponding to an intrinsic luminosity near 1 x 10 40 ergs _1 establishing it as an ultraluminous X-ray 
source (ULX). The energy spectrum is dominated by disk emission with a weak or absent Compton 
component and there is no significant timing noise; both are indicative of the thermal state of stellar- 
mass black hole X-ray binaries and inconsistent with the Compton-dominated state typical of most 
ULX spectra. A previous measurement of the X-ray spectrum shows a harder spectrum that is well 
described by a powerlaw. Thus, the binary appears to exhibit spectral states similar to those observed 
from stellar-mass black hole binaries. If the hard state occurs in the range of luminosities found for 
the hard state in stellar-mass black hole binaries, then the black hole mass must be at least 85Af Q . 
Spectral fitting of the thermal state shows that disk luminosities for which thin disk models are 
expected to be valid are produced only for relatively high disk inclinations, > 60°, and rapid black 
hole spins. We find a* > 0.98 and M > 154 M Q for a disk inclination of 60°. Higher inclinations 
produce higher masses and somewhat lower spins. 

Subject headings: black hole physics — galaxies: individual (blue compact dwarf, I Zw 18) — X-rays: 
binaries — X-rays: galaxies 



1. INTRODUCTION 

The early universe was deficient in heavy elements. 
Thus, understanding the metallicity dependence of X- 
ray emission and compact object formation is impor- 
tant for our understanding of the earl y universe, particu- 
larly the mechanism s of reionization (|Madau et al.l 120041 : 
iMirabel et al.ll20TlT) and the effect of feedback from com- 
pact objects prod uced by star forma tion on the dynamics 
of early galaxies (jPower et al.ll2009l ). 

Recent models of the formation and evolution of stars 
and binaries with low metallicity suggest that their 
X-ray emission is significantly enhanc ed (|Dravl 120061 : 
iDrav fe Tout] 120071 iLinden et all |201(I. These trends 
are supported by observations. iMapelli et al.l (|2010f ) 
demonstrated an increase in the number of X-ray bright 
HMXBs ( for a given SFR) for moderately low metallicity 
galaxies. IKaaret et al.l (|2011[ ) found a suggestion of en- 
hanced X-ray luminosity (for a given SFR) in a sample 
of very-metal-poor blue compact dwarf galaxies. 

It has also been suggested that low metallicity envi- 
ronments may lead to the formation of unusually mas- 
sive compact objects. Solar metallicity stars are ex- 
pected to lose much of their mass via winds at the 
ends of their lives and produce remnants with masses 
less than about 20 Mq. Reduced metallicity decreases 
the effect of winds and leads to more massive remnants 
(jBelczvnski et al.ll2010f ). This has been suggested to ex- 
plain the increased occurrence of X-ray binaries with 
very high luminosities, the ultralumino us X-ray sources 
(|Kaaret et al.l l2001: Fen g fc Soriall2011ft. in low metallic- 
ity environments (jPakull fc Mirionill20 03: Zampicr i et al.l 
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I200J IMapelli et al.ll200l . 

Blue compact dwarf galaxies (BCDs) are physically 
small g alaxies with blue opt i cal colors and low metal- 
licities (jKunth fc Ostlidl2000l : IWu et al.ll2006h . They are 
the best local analogs to early galaxies yet identified. I 
Zw 18 is a nearby BCD with the lowest known metallic- 
ity and the highest X-ray luminosity (Kaaret et al. 2011). 
The X-ray emission from I Zw 18 is dominated by a sin- 
gle X-ray binary and observations reported to date have 
shown a maximum luminosity below the Eddington lu- 
minosity of a 20 Mq compact object. 

Here, wc report on a measurement of the X-ray spec- 
trum of I Zw 18 that reveals the highest flux yet observed. 
The spectrum shows distinct curvature indicative of a 
black hole X-ray binary in the thermal state. We present 
our results on the X-ray spectrum in §2 and discuss their 
implications for the nature of compact object in §3. 

2. OBSERVATIONS AND ANALYSIS 

I Zw 18 was observed by the X-ray Multi-Mirror Mis- 
sion (XMM-Newton) on 2002 April 10 for 32.8 ks and on 
2002 April 16 for 28.9 ks. We reduced the data according 
to the standard procedures for imaging spectroscopy with 
the European Photon Imaging Camera (EPIC). Back- 
ground flaring reduced the useful exposure for the first 
observation for the pn camera to 24.1 ks and for the two 
MOS cameras to 30.9 ks. Flaring was much worse for 
the second exposure and reduced the useful exposure for 
the pn to 5.3 ks. Thus, we chose to further analyze only 
the first observation. 

We extracted spectra using circular source regions with 
30" radius centered on the X-ray source apparent in each 
image. This extraction region is larger than the extent 
of I Zw 18 and therefore includes all X-ray emission from 
the galaxy. Background subtraction was performed using 
a circular background region with a radius of 60" located 
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Table 1 

X-ray Spectral Fits 



Model 



X 2 /D°F 



Lx 

(10 40 erg s- 1 ) 



(10 2 



N H 



kT/E c 
(keV) 



Powerlaw 

Diskbb 

Cutoff powerlaw 

Diskbb+powerlaw 

Simpl*kcrrbb 



220.8/151 
172.0/151 
156.7/150 
157.3/150 
153.4/150 



3.0 
1.1 

1.2 
1.3 
1.2 



3.8±0.5 




2.01±0.06 


0.7±0.3 


1.05±0.05 




1.4±0.6 


2.1±0.6 


0.80±0.26 


1.4±0.6 


1.04±0.07 


2 


1.2±0.4 




2 



Note. — The table includes: the m odel name, goodn ess of fit (x ) and degrees of freedom (DoF), the intrinsic luminosity in the 0.3—10 keV 
band assuming a distance of 18.2 Mpc (|Aloisi et al.l|2007|) and isotropic emission, absorption column density for the component with abundances 
fixed to Z/Zq — 0.019, the disk temperature at the inner edge (kT) or the powerlaw cutoff energy (E c ), and the photon index (T). 
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Figure 1. X-ray spectrum of I Zw 18. The solid, stepped curves 
show the best fit to the simpl*kcrrbb model with absorption as 
described in the text. There are three curves and three sets of 
data points. The upper curve/data is for the pn. The lower two 
curves/data are for the MOS1 and MOS2. 



on the same CCD chip. Together the 3 spectra contain 
about 3940 net counts. Spectral response files suitable 
for point source analysis were calculated using the most 
recent calibrations. The spectra were grouped to have 
a minimum of 16 counts per bin. We fi tted the X-ra y 
spectra using the xspec software package (|Arnaudlll996[ ). 
We used the 0.3-10 keV energy range for fitting. 

The Chandra X-ray Observatory observed I Zw 18 for 
41 ks on 2000 February 8 an d imaging and s p ectro - 
scop ic results were repo rted by iBomans fc Weisl (|2002[ ) 
and iThuan et al.l (|2004D . The imaging shows that the 
X-ray emission was dominated by a single point source 
with at least 96% of the flux. Both groups found good 
fits with an absorbed power-law model, but residuals 
near 0.65 keV that th e y ascr ibed to an OIII hydrogen- 
like line. IThuan et al.l (|2004| ) report an observed flux of 
7.2 x 10 -14 erg cm~ 2 s _1 in the 0.5-10 keV band, a photon 
index of2.01±0.16, and an absorption column density, 
N H = (1.44 ± 0.38) x 10 21 cm- 2 . 

We fitted the XMM spectra using an absorbed power- 



law model and found an observed flux of (2.78 ± 0.12) x 
10~ 13 ergcm~ 2 s _1 in the 0.5-10 keV band, a photon in- 
dex of 2.31 ± 0.09, and an absorption column density, 
N H = (2.7 ±0.3) x 10 21 cm -2 . Allowing the normaliza- 
tion to vary between the 3 detectors did not improve the 
fit, so we used the same normalization for all three detec- 
tors during all subsequent fitting. The marked increase in 
flux from the Chandra to the XMM observation shows 
that the emission is variable and likely predominantly 
due to a single X-ray binary. Thus, we consider only 
spectral models appropriate for X-ray binaries. Since 
the XMM extraction region includes the whole galaxy, 
it is possible that a second source, other than the one 
detected with Chandra, contributes a significant fraction 
of the flux observed in the XMM observation. Due to the 
low numbers of X-ra y sources detected in similar galaxies 
(jKaaret et al.ll20TTT) . this is unlikely, but could be tested 
with a new Chandra observation. 

The absorption column density required for these spec- 
tral fits is well above the total Galactic Hi column den- 
sity towards I Zw 18 of N H = 2.5 x 10 20 cm" 2 . Thus, 
most of the absorbing material likely resides in I Zw 
18. The metallicity of I Zw 18 is measured via op- 
tical spectroscopy of H n regions to be Z/Zq = 0.019 
(|Izotov fc Thuanl 1 19990 . where we have adopted a so- 
lar oxygen abundance of 12+log(0/H) = 8.9, while the 
abundances in th e neutral interstel lar medium are sev- 
eral times lower (jAloisi et al.l 120031 ). The X-ray binary 
imaged with Chandra lies in a star formation region and 
the Hn region abundance is appropriat e. The X-ray ab- 
sorpti on m odel used above an d also by IBomans fc Weisl 
(|2002D and IThuan et al.l (|2004D assumes solar abundance 
and is likely an incorrect description of the true energy 
dependence of the absorption. In the fits below, we use 
two absorption components: one with solar abundances 
with a column density fixed to the Galactic Hi column 
towards I Zw 18 (TBabs in xspec) and a second with 
a variable column density (TBvarabs) , abundances fixed 
to Z/Zq = 0.019, and redshift fixed to 0.00254. 

The results of fitting to various models are shown in 
Table Q] The only model which is excluded is the simple 
powerlaw, while the diskbb model is marginally excluded. 
A model consisting of the sum of these two components 
provides a good fit. The three best fitting models are sta- 
tistically indistinguishable. They have an observed flux 
of 2.7 x 10~ 13 crgcnr 2 s" 1 in the 0.3-10 keV band and 
the absorption column de nsity within I Zw 18 is in the 
range 1.2-1.4xl0 21 cm- 2 . ILeili et al.l ([2012D mapped Hi 
in I Zw 18 via the 21 cm line at 2" resolution and their 
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Figure 2. Dependence of black hole mass (Mbh, upper panel) 
and Eddington ratio (L/L^dd, lower panel) on black hole spin (a*) 
for disk inclinations of i = 60° (circles) and i = 75° (diamonds). 
The dashed line in the lower panel shows L/L^dd = 0.3. 

maps show a total column density near 6 x 10 21 cm~ 2 
at the X-ray source position. The lower Njj values from 
the X-ray spectral fits are reasonable if the X-ray bi- 
nary lies closer than the midplane of the galaxy or if the 
metallicity of the Hi gas is lower than that of the Hn re- 
gions, suppressing X-ray absorption. We note that none 
of these model fits requires an emission line at 0.65 keV 
and we suggest that the line was due to an inappropriate 
choice of absorption model. 

The cutoff powcrlaw model is empirical, but provides 
a simple analytical form, f(E) = E~ r cxp(— E/E c ). The 
improvement in the fit going from the powerlaw to the 
cutoff powerlaw indicates spectral curvature at high en- 
ergies. 

The sum of a powerlaw and multicolor disk blackbody 
(diskbb) is often used to mo del the spectral of Galac- 
tic b lack hole X-ray binaries (|Remillard fe McClintockl 
120061) . In our fitting, the photon index (T) was not well 
constrained, so we fixed T = 2. The powerlaw normal- 
ization was 2.3 x 10 -5 and the diskbb normalization was 
7.8 x 10 -3 . The disk component produces 76% of the 
flux in the 2-10 keV band. 

The kcrrbb model represents a sophisticated model of a 
thin accretion disk around a Kerr black hole including all 
relativistic effects and also self-ir radiation of the disk due 
to light deflection (|Li et al.H2005D . The simpl model adds 
a Comptonization co mponent calculated by convolution 
of the disk spectrum dSteiner et al.ll2009D. W e fixed the 
source distance to 18.2 Mpc (jAloisi et al.ll2007f). the spec - 
tral hardening factor to 1.7 dShimura fe T akahara 19951 ). 
the torque at the inner disk boundary to zero, and in- 
cluded self irradiation and limb darkening in the model. 
The fit was insensitive to the black hole spin param- 
eter (a*), inclination (i), and Comptonization photon 
index (r). The best fitted spectrum for a* = 0.9986, 
i = 60°, T = 2.0 is shown in Fig. [Q The fitted 
parameters not included in Table [1] were: black hole 



tion, 0.27 ± 0.13. We also performed fits with T = 2.6, 
to cover the range typically found in the thermal state, 
and none of the fit parameters changed significantly. 

We performed fits for spins over a range from a* = 0.9 
to th e maximum achievable spin a* = 0.9986 (|Li et al.l 
120051 ) for inclinations of 60° and 75°. Lower inclinations 
and lower spins produce Eddington ratios, L/L^dd > 0.3, 
inconsistent with application of the kcrrbb model. The 
spin was kept fixed for each individual fit. The x 2 var- 
ied by less than 0.5 over this range. The Comptonized 
fraction did not vary significantly, with best fitted val- 
ues in the range 0.26-0.33. The absorption column also 
did not vary significantly, with best fitted values in the 
range 1.0-1.2 xlO 21 cm -2 . The best fitted black hole 
mass (Mbh) and accretion rate vary with both a* and i. 
Fig. [2] show the variation in Mbh and the ratio of disk 
luminosity to the Eddi ngton luminosity {L/ L p,dd)- 

For comparison with iGladstone et al.1 (|2009f ) , we fitted 
the XMM data in the 2-10 keV band with a powerlaw 
and a broken powcrlaw, each without absorption. We 
note that our spectrum cont ains fewer counts , ~ 40 00, 
than the > 10, 000 required by IGladstone et al.1 (|2009D for 
inclusion in their sample. The powerlaw provides an ad- 
equate fit with x 2 /DoF = 51.7/43 and a photon index 
T = 2.40 ± 0.16. The broken powerlaw provides only 
a slight improvement, x 2 /DoF = 48.4/41 for an F-test 
value of 0.25. The fit is insensitive to the low energy pho- 
ton index, so we fixed it to 1.4 (one less than T for the 
simple powerlaw). The break energy is then between the 
lower end of the fitting range, 2.0 keV, and 2.9 keV and 
the high energy photon index is 2.59±0.25. We note that 
fitting the full energy range to a broken powerlaw model 
with absorption leads to a break energy of 1.9 ± 0.3 keV. 
We also fitted the data in the full energy range with 
the sum of a disk model (diskpn) and a Comptonization 
model (comptt) with absorption. We tied the Comp- 
tonization photon input temperature to the disk temper- 
ature and used an inner disk radius of QRq. The model 
produced a reasonable fit with x 2 /DoF = 152.6/148 and 
a disk temperature kT = 0.51q2 keV. However, the other 
parameters were very poorly constrained. Only lower 
bounds were obtained for the Compton optical depth, 
t > 0.05, and th e plasm a temperature, kT e > 1.14keV. 
IGladstone et al.l ([2009) found significant improvement 
(F-test significance level of >99%) of the best fitted 
model relative to a hot corona with kT e = 50 keV for the 
spectra of sources identified as in the ultraluminous state. 
Making the same comparison, we find A% 2 = 0.8 (DoF 
= 148) and an F-test significance level of 62%. Thus, 
the low temperature Compton component does not pro- 
duce a statistically significant improvement in the fit, as 
required for the ultraluminous state. 

We reduced the Chand ra d ata previously ana lyzed 
by iBomans fe Weii ([2003 ) and IThuan et all ([200l and 
found a net 490 counts. Fitting with an absorbed pow- 
erlaw, we found resu l ts co nsistent within errors with 
those of IThuan et al.l (|2004f ) . For comparison with the 
XMM-Newton results reported above, we also fitted the 
Chandra data with models with two absorption compo- 
nents, a Galactic component with solar abundances and 
Nh = 2.5 x 10 20 cm -2 and a component for I Zw 18 
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Figure 3. X-ray spectra of I Zw 18 at low and high flux levels. The 
points and the solid curve show the XMM observation fitted with 
the cutoff powcrlaw model. The dashed curve shows the powcrlaw 
model fitted to the Chandra data. 

with Z/Zq = 0.019 with variable Nh- Fitting with a 
powerlaw, we found a good fit with % 2 /DoF = 16.2/25, 
r = 1.78 ± 0.21, N H = §tf x 10 20 cm- 2 , and a flux of 
8.3 x 10- 14 ergcm- 2 s" 1 in the 0.3-10 kcV band. Fitting 
with a cutoff powerlaw produced an essentially identi- 
cal fit with a best fitted cutoff energy of 500 keV and 
a 90% confidence lower bound on the cutoff energy of 
4.7 keV. Fitting with a diskbb model produced a poor 
fit with x 2 /DoF = 36.7/25. Thus, we conclude there is 
no evidence for spectral curvature at high energies in the 
Chandra spectrum. 

We examined the timing properties of the XMM pn 
data. The 0.3-10 keV light curve binned in 450.8 s in- 
tervals appears constant, with x 2 /DoF = 55.8/54. The 
rms power integrated over 0.01-1 Hz is less than 5% at 
99% confidence. Due to the low number of counts, the 
Chandra data do not allow us to place useful constraints 
on the timing properties. 

3. DISCUSSION 

The XMM-Newton observations described here show 
the highest flux yet observed from I Zw 18, correspond- 
ing to an intrinsic luminosity near 1 x 10 40 ergs _1 and 
establishing the X-ray binary in I Zw 18 as an ultralu- 
minous X-ray source (ULX). 

The spectral shapes of the XMM versus Chandra data 
arc different with the high flux (XMM) spectrum show- 
ing distinct curvature at high energies while the low flux 
(Chandra) spectrum shows no evidence of curvature, see 
Fig. [3] This is consistent with the state transitions seen 
in stellar-mass black hole X-ray binaries, specifically the 
transition between the hard X-ray spectral state and the 
thermal state. The lack of timing noise and low frac- 
tion (< 30%) of powerlaw flux in the diskbb plus pow- 
erlaw model reported in Table [1] for the high flux state 
are consistent with its id entification as the thermal state 
(iRemillard fc McClintock1[2006l ). 

iGladstone et al.l (|2009f) has suggested the existence of 
an 'ultraluminous' state based on curvature at high ener- 
gies in the spectra of several ULXs. Broken powerlaw fits 
to the XMM spectra of I Zw 18 produce break energies 



well below any of those reported for the ultraluminous 
state. Also, a spectral model with a low temperature 
Compton component {kT e = 1-3 keV) produces no sta- 
tistically significant improvement in the fit relative to a 
model with a hot Compton component {kT e =50 keV). 
These results suggest that the XMM spectrum is incon- 
sistent with the criteria established for the ultraluminous 
state. However, we caution that the spectrum contains 
fewer counts than those used to establish the properties 
of the ultraluminous state. A higher quality spectrum 
with at least 3 x the number of counts would be needed 
to draw definitive conclusions. 

State transitions have been reported in ULXs that 
have been modeled and interpreted as being un- 
like those of stellar-mass black hole X-ray binaries. 
(jPintore fe Zampieri|[20T2T) analyzed multiple XMM ob- 
servations of the two ULXs in NGC 1313 and classified 
their spectra into two states: the 'thick corona' state and 
the 'very thick corona' state. The spectra in both states 
are dominated by a coronae with kT e in the range 1— 
6 keV. Other authors have reported on spectral evolution 
in ULXs with the common characteristic that a Comp- 
ton component with kT e >1 kcV, cont ributes a major 
and usually dominant fra ction of the flux dFeng fe Kaaretl 
I2009t iKaiava et aT1l2012f ). In the XMM spectra of I Zw 
18, a Compton component modeled as a powerlaw pro- 
duces a low fraction of the X-ray flux and a Compton 
component with a low temperature (kT e = 1-3 keV) pro- 
duces no statistically significant improvement in the fit. 
Thus, the spectral state transitions in most ULXs appear 
different from that seen in I Zw 18. 

While rare, the thermal state and transitions between 
the hard state and the thermal state are sometimes found 
in ULXs. The strongest case for such a transition is for 
M82 X-l which shows both a spectral transition and a 
simultaneous change in the X-ray timing properties of 
the so urce as seen in stellar -mass black hole X-ray bi- 
naries (|Feng fc Kaaretl l2010f) . Spectral evidence for the 
hard/thermal state transition has also been presented 
for the extremely lum inous X-ray binary in E SQ 243-49 
(IServillat et al.ll2C>Tlh and for M82 X37.8+54 (IJin et al.l 
I2010D . Spectral e vidence for the t hermal state was found 
for NGC 247 X-l (|Jin et al.H20Tlh and a t lower luminosi- 
ties in CXOM31 J004253. 1+411422 (jMiddleton et al.l 

HoU). 

Detection of a state transition similar to those seen 
in stellar-mass black hole X-ray binaries strengthens in- 
terpretation of the hard spectrum seen in the Chan- 
dra observation as evidence that the source was in 
the hard state. The maximum luminosities observed 
from stellar-mass black holes in the hard state are 
less than 0.3Lprir\ ([Rodriguez, Corbel, fc Tom sick 2003; 
Zdziarski et~aTll2004t lYuan et al.l I2007t IMivakawa et all 



20081 ). The luminosity of the binary in I Zw 18 while in 
the hard state was 3.3 x 10 39 ergs _1 . This corresponds 
to L < 0.3L E dd only if M BH > 85M . 

In the thermal state, the inner radius of the accre- 
tion disk is set by the mass and spin of the black hole. 
As described above, we found that the spectral fits with 
simpl*kerrbb model were insensitive to the black hole 
spin (a*) and inclination (i), thus, we considered a range 
of values for a* and i. For i < 60°, the source is 
supcr-Eddington for all spins. The accretion disk is ex- 
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pected to be geometrically thin, and thus the kerrbb 
model is expected to b e valid, only for L/LEdd < 0.3 
(jMcClintock et aLll2006f ). so we considered only fits with 
i > 60°. Fig. [2] shows the dependence of black hole 
mass (Mbh) and disk luminosity relative to Eddington 
(L/L Edd ) on a* for i = 60°. We find that L/L Edd < 0.3 
(within errors) for a* > 0.98 with M BH = 195i"M . 
For higher spins, L/LEdd decreases and Mbh increases. 
For i = 75°, we find that L/L Edd < 0.3 for a* > 0.92 
with Mbh = 288^57-M©. Again, L/L E dd decreases and 
Mbh increases for higher spins. 

In conclusion, the X-ray spectrum from the binary in 
I Zw 18 while in the high flux state can be interpreted in 
terms of a non-rotating black hole with super-Eddington 
emission and a mass similar to that of the most massive 
known stellar-mass black h ole (|Silverman fc Filippenkq 
120081: iPrestwichetH] 120071) . However, if the nature of 
the state transition seen from the source is similar to 
those seen from stellar-mass black hole X-ray binaries, in- 
cluding the luminosity threshold for the transition, then 
the compact object is likely a near maximally-rotating 
black hole with an unusually high mass Mbh > 85Af Q . 
This is close to the maximum mass of 75M predicted for 
black holes formed via stellar col lapse in a low metallicit y 
(Z/Zq = 0.019) environment (jBelczvnski et alJl2010T >. 
Modeling of the thermal state spectrum suggests a higher 
mass, M B h >= 154M Q , in the IMBH range. We note 
that the X-ray source is coincident with a massive star 
cluster, sites suggested as potential places of original for 
IMBHs (jPortegies Zwart et"aLll2004D. particula rly in low- 
mctallicity environments (jMapelli et al.1 120131 ). Due to 
their similarities to early galaxies, identification of an 
IMBH in a BCD would be of particular interest for our 
understanding of the fo rmation of early ga laxies and su- 
permassive black holes (|Reines et al.ll201ll ). 
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